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Formulation of a unified method for low- and high-energy 
expansions in the analysis of reflection coefficients 
for one-dimensional Schrodinger equation 

Toru Miyazawa®') 

Department of Physics, Gakushuin University, Tokyo 171-8588, Japan 

We study low-energy expansion and high-energy expansion of reflection coefficients 
for one-dimensional Schrodinger equation, from which expansions of the Green func¬ 
tion can be obtained. Making use of the equivalent Fokker-Planck equation, we 
develop a generalized formulation of a method for deriving these expansions in a uni¬ 
fied manner. In this formalism, the underlying algebraic structure of the problem can 
be clearly understood, and the basic formulas necessary for the expansions can be 
derived in a natural way. We also examine the validity of the expansions for various 
asymptotic behaviors of the potential at spatial infinity. 


I. INTRODUCTION 


Our object of study in this paper is the one-dimensional steady-state Schrodinger equation 


+ Vs{x)'4){x) = k'^ipix), 


( 1 . 1 ) 


where Fs is a real-valued function, and fc is a complex number in the closed upper half plane 
(Im/c > 0). We assume that the origin of the energy scale = 0) is set at the energy of the 
ground state. [When dealing with a Schrodinger operator —(P/dx^-\-Ua,{x) which has a ground state 
with energy Eq ^ 0, the ground-state energy can be shifted to zero by defining Vs (a;) = Us{x) — Eq] 
In other words, corresponds to the energy relative to the ground state. The ground state may 
either be a bound state or a continuum state. 

It is well known that (1.1) is equivalent to the Fokker-Planck equation^ 


- dl 




(l){x) \ = k‘^(j){x), 


( 1 . 2 ) 


which describes diffusion in a potential V(x). We define 

/(a:) =(1-3) 

Equations (1.1) and (1.2) are related by 

^(x) = Vs{x) = [f{x)]'^ -b ^/(^)- (1-4) 

Substituting (j) = into (1.2) yields Schrodinger equation (1.1). It is also easy to see that 

tpo{x) = is the ground-state wavefunction satisfying (1.1) with k = 0. 

Let U{x,y;k) be the 2x2 matrix satisfying the differential equation 




f-ik f{x)'\ 
\f{x) ik J 


U{x,y,k) 


(1.5) 


with the boundary condition U(y,y;k) = I (identity matrix). Here, / is the the function defined 
by (1.3). The elements of U can be written with two functions a and f3 as 


U{x,y;k) 


a{x,y;k) I3{x,y;-k)\ 
I3{x,y;k) a{x,y,-k) J 


( 1 . 6 ) 
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We define the transmission coefficient r, the left reflection coefficient i?/, and the right reflection 
coefficient Rr as 


T{x,y;k) 


1 

a[x,y\k)' 


Ri{x,y, k) = 


^{x,y]-k) 

a{x,y;k) 


Rr{x,y;k) 


I3{x,y;k) 
a{x,y; k)' 


(1.7) 


(The reason why they are called the transmission and reflection coefficients is explained in Ap¬ 
pendix A of Ref. 2.) If V{x) is sufficiently well behaved at infinity, we can define the reflection 
coefficients for semi-infinite intervals as 


Rr{b,—co;k)= lim Rr{b,a;k), Ri{co,a;k)= lim Ri{b,a;k). (1-8) 

a—> —oo b—J--I-00 

When fc is a real number, it may happen that the limits in (1.8) do not exist. In such cases, we 
interpret (1.8) as the limit approached from the upper half plane, 

i?r(^, —oo; A) = lim lim Rr{b,a;k + ie), i?;(oo, a; fc) = lim lim Ri{b,a;k + ie). (1-9) 

e4-0 a^ — oo €4-0 b—>--\-oo 


Let Gs{x,y;k) be the Green function of Schrodinger equation (1.1) satisfying 




- Vs (a;) + fc 


Gsix,y;k) = S(x - y) 


( 1 . 10 ) 


with the boundary conditions Gs(x,y;k) —>■ 0 as jx — yj —>■ oo for Imfc > 0. [For Imfc = 0, we 
define Gs(a;, y; k) = lirngj^o Gs(a;, y,k -\- ie)] We can express this Gs as^ 


Gs{x,y;k) = — 


1 


: exp 


r 

ik{x — y) — ik / S{z,k)dz 
Jv 


where we have assumed (without loss of generality) that x > y and defined 

S{x, k) = Sr{x, k) + Si{x, k), 


( 1 . 11 ) 


( 1 . 12 ) 


Srix, k) 


Rr{x, —oo; k) 

1 + Rr(x, —oo; k) ’ 


Si{x,k) 


Ri{co; x; k) 

1 -I- Ri{co, x; k) 


(1.13) 


Asymptotic behavior of solutions of the Schrodinger equation in low- and high-energy regions has 
been an important object of study for many years.Since the Green function can be expressed 
solely in terms of Rr{x, —oo; k) and Ri(oo, x; k) as shown in (1.11), all the information about 
the behavior of solutions can be obtained through the analysis of the reflection coefficients. The 
study of the reflection coefficients is also important in spectral theory and inverse problems. In 
particular, Sr{x,k) and Si{x,k) defined by (1.13) are closely related to the Weyl-Titchmarsh m- 
functions, which play a significant role in spectral analysis of the Schrodinger operator.(The 
relation between Sr, Si and the m-functions is discussed in Appendix C of Ref. 2.) In this paper, 
we investigate low- and high-energy expansions of Rr{x, —oo; k). The results for Ri{co,x;k) can 
be obtained in the same way. 

To study low- and high-energy expansions, the Fokker-Planck equation is more suitable than 
the Schrodinger equation itself because the structure of the problem becomes more transparent 
for the Fokker-Planck equation. By using the Fokker-Planck potential V, we can carry out the 
analysis more systematically, and more explicit expressions can be obtained, than by using Vs- In 
this paper, we deal with expressions in terms of V{x) or f{x). 

We assume that V{x) is a piecewise continuously differentiable^^ real-valued function. Our 
method is applicable to V(x) with various asymptotic behaviors at spatial infinity, including the 
cases where V{x) is finite or ±oo as a; —>■ ±oo. Gorresponding to such V^(a;), the Schrodinger 
potential Vs{x) is finite or -|-oo as a; —>■ ±oo. [We do not consider the cases where lima;^±oo Vs(a^) = 
—oo.] We will also deal with asymptotically periodic potentials. Specific conditions on the potential 
will be given when they become necessary. 

In the formalism based on the Fokker-Planck equation, we can express Rr{x, —oo; k) in terms 
of a linear operator. The high- and low-energy expansions of the reflection coefficient then reduce 
to the expansions of this operator. This method has been studied in the previous papers. 
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However, the derivation of the expansions in these papers was not rigorous. The origin of the 
basic formulas used in these works was not clear, and the meaning of the method has not been 
fully understood. It is the aim of the present work to construct a generalized formulation which 
provides a clear overall view of the problem and which gives a rigorous justification to the intuitive 
arguments of the previous papers. In this generalized framework, we will clarify the mathematical 
structure of the method and show how these expansions, together with the intermediate formulas, 
can be derived in a simple and natural way. 


II. PRELIMINARIES 


We can write Eq. (1.5) in the form 


d 


—U (x, y] k) = [2ikJ3 - 2f{x)Ji]U{x, y; k) 


with 




1 /O I 


2(10 


J 2 = 


1 /O -i 


2 Vi 0 


J3 = 


1 f-l 0 
0 1 


The matrices Ji, J 2 , J 3 satisfy the commutation relations 

[</2, ds] = *<^ 1 , [d3,Jl]=iJ2. 

We also define J+ = Ji + i J 2 and J_ = Ji — i J 2 . Then, U(x, y; k) can be expressed as 
U{x,y;k) = exp[-Rr{x,y; k)J+] [r(x, j/; exp[Ri{x,y; k)J-], 


( 2 . 1 ) 


( 2 . 2 ) 


(2.3) 


(2.4) 


where [r(x, y; = exp{[2 logT(a;, y; A:)] J 3 }. Indeed, by using explicit forms (2.2), we can easily 

verify (2.4) by directly calculating the right-hand side as^® 




1 0 
Rr 1 


T-l 0 

0 T 


1 -Ri\ _ ( a{k) (i{-k) 

0 1 y V 


(2.5) 


Actually, we can prove (2.4) only by using commutation relations (2.3), without reference to explicit 
expressions (2.2) of Ji. This means that (2.4) is a representation-independent expression.^® Let 
Ji, J 2 , J 3 be any set of operators satisfying commutation relations (2.3). Then, the solution of 
Eq. (2.1) with the boundary condition U{y,y) = 1 (identity operator) is given by (2.4). 

Treating / as a perturbation, we can derive from (1.5)-(1.7) the expressions of r, Ri, and Rr as 
formal series in powers of /, 


_ ^ik(x-y) _ ^ik(x-y) 


r dz, r dz 2 fizi)fiz 2 ) +. 

'Jy 'Jy 


^2ik{zi-y) 


T{x,y;k) = e* 

Ri{x,y]k) = - dzif{zi)< 

Jy 

- r dz^ r dz2 r dzsf{z,)f{z2)f{zs)e^^'^^^^-^^+^^-y^ 

J V j V j Zo 


> Z2 
^2ik{x-zi) 


Rr{x,y]k)= / dzif{zi)e 
Jy 

+ r dz^ r dZ2 r dZ3f{z,)f{z2)f{z3)e^^’^^--^^+^^-^^ 
J y J Z^ J V 


( 2 . 6 a) 


( 2 . 6 b) 


( 2 . 6 c) 


Equations (2.6) can be represented by diagrams as in Fig. 1(a). These diagrams are to be in¬ 
terpreted according to the rules shown in Fig. 1(b). Each line segment connecting the points xi 
and X 2 corresponds to exp[ifc(x 2 — xi)], and a factor ±f{z) is assigned to each turning point z. 
Integration over the positions of Zi is implied in Fig. 1(a). (The vertical direction of the diagram 
does not have any meaning.) 

Such diagrams are a useful tool for understanding the structure of the transmission and reflection 
coefficients in an intuitive way. The series in (2.6) are formal, and we need not be concerned about 
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(a) 


ib) 




FIG. 1. (a) Diagrammatic representation of Eqs. (2.6). (b) Rules for interpreting the diagrams. 


their convergence. Although we shall make use of diagrams in Sec. Ill, the obtained results [such 
as (3.5) and (3.6)] are valid even if / is not small and the series in (2.6) are not convergent. These 
results can always be proved directly by using (1.7), without using series expressions (2.6). 

The transmission and reflection coefficients become easier to handle by introducing an addi¬ 
tional variable as follows.^® Suppose that V{x) has a jump of magnitude w &t x = xq, i.e., 
lime 4 ,o [I^(a^o + e) ~ V(xo — e)] = w. Then, f(x) has a delta function singularity as 

fix) =-^S(x - xo)-{ -. (2.7) 

From (1.5) and (2.7), we find that the matrix U across the discontinuity takes the form 


lim U (xo + e, xq 


e: k) = lim exp 
^ €4.0 




dx 



—wf2 

0 


fS{x - Xo) 

( cosh f 
" sinh f 


-^5{x - Xo) 
ik 

-sinhf\ 

cosh^ ) 


( 2 . 8 ) 


This matrix, which does not depend on xo, describes the jump of the Fokker-Planck potential V. 
We multiply U{x,y; k) [Eq. (1.6)] from the left by this matrix and define 


f a{x,y,w,k) I3{x,y;w;-k)\ ^ f cosh f - sinh f \ . 

y^{x^y;w;k) a{x,y;w;—k) J \^ — sinh ^ cosh ^ J v ’y’ / 


(2.9) 


(The bar does not mean complex conjugation.) Using a and (3 we define, as we did in (1.7), 
1 


r = 


a{x,y;w;k)' 


l3{x,y;w;-k) - _ I3{x,y;w;k) 

— _ / , \ 1 -^r — 


a{x,y;w;k) ’ 


a(x, y; w; k) 


( 2 . 10 ) 


We can interpret f, Ri, and Rr as the transmission and reflection coefficients that include the effect 
of a jump in V(x) at the right endpoint of the interval. The additional variable w corresponds to 
the magnitude of the jump. It is convenient to introduce, instead of w, the new variable 


^ = tanh 


w 


( 2 . 11 ) 
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and regard r, Rr, and Ri as functions of {x,y,^,k}. For simplicity, we shall often omit to write 
the argument k. From (1.6), (1.7), (2.9), and (2.10) we have^® 


r{x,y;0 


y^l-^^T{x,y) 

l-^Rr{x,y) 


Rr{x,y;0 


Ri{x,y]0 = Riix,y) + , n 

1 - iRr[x,y) 

(1 - ^^)Rr{x,y) _ 

l-^Rr{x,y) 


( 2 . 12 ) 


III. STRUCTURE OF THE GENERALIZED TRANSMISSION AND REFLECTION COEFFICIENTS 


Let us further generalize (2.12) and define 

yT{x,y) 


T{x,y;^,fj,) = 


l-^Rrix,y)’ 


Lix,y;0 = Ri{x,y) + 


ix'^{x,y) 

1 - ^Rr{x,y)'' 


f,, c f^'^Rrix,y) 

R[x,y;^,y) = 


1 - ^Rr{x,y)’ 

where ^ and y are complex variables. Equations (2.12) can be written as 

x{x,y,^) = f{x,y,^,'y), Ri{x,y;^) = L{x,y;^), 

Rr{x, y; 0 = Rix, y; C, 7 ) - C, 


(3.1) 


(3.2) 


where 7 = ■ While ^ and 7 are related by + 7^ = 1, the variables ^ and y in (3.1) are 

independent. Moreover, we regard ^ and y as complex numbers. The original r, Ri, and Rr are 
recovered from T, L, and R, respectively, by setting ^ = 0 and y = 1. 

The reflection coefficient has the property that |i?r| < 1 for Im/c > 0. So, obviously T, L, and R 
are analytic functions of ^ in the unit circle |^| < 1. The meaning of Eqs. (3.1) can be intuitively 
understood by expressing them as infinite series 

00 00 00 

f = yTY,mrr, L = Ri+iT^Y.^iRrT, R = y'^RrY,mrT, (3.3) 

n—0 n—0 n—0 

which are convergent for |^| < 1. These equations are illustrated in Fig. 2. It can be seen from the 
definition that t{x,x) = 1, Rr{x,x) = Ri{x,x) = 0. Hence, 

f(x, x; C, y) = y, L{x, x; ^ R{x, x; y) = 0. (3.4) 

Let z < y < X. In the diagrammatic representation, T(x, z), L{x, z), R(x, z) can be constructed 
from T(x,y), L{x,y), R(x,y) and T{y,z), Ri{y,z), Rr{y,z) as shown in Fig. 3. Comparing Fig. 3 
with Fig. 2, we find that T{x, z) and L(x, z) are obtained from T(y, z] y) and L(y, z] by putting 
L(x,y) and T(x,y) in place of ^ and y, respectively. That is, 

T{x, z; y) = f{y, z; L{x, y; O^T^x, y; y)), L{x, z; ^ = L{y, z; L{x, y; ^)). (3.5) 


In the bottom row of Fig. 3, the first term on the right-hand side is exceptional and must be treated 
separately. So, the expression for R{x, z) includes an additional term as 

R{x, z; C, y) = R{y, z; L{x, y; ^), f{x, y; y)) + R{x, y; y). (3.6) 

Differentiating (3.5) and (3.6) with respect to x, and then setting y = x, we obtain 


9 •' 9 - 9 - 

—T{x,z;^,y) = —T{x,z-,^,y) —L{x,y]^) 

^L(x,z;C) = ^L(x,z;C) -^L{x,y;^) 


y^x 


9 - 9 - 


y^x 


y^x 
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FIG. 2. Schematic illustration of Eqs. (3.3). 


z; fj.) = -^R{x, z] C, m) tT^^x, y; C) + -^Rix, z; y) t-T^x, y, y,) 
dx dS, dx dy dx y^^ 

+ -^R{x,y,i,y-) , (3.7) 

y —21 

where we have used (3.4). From Fig. 4, we find that (see Appendix A) 
d - 

-7r'Rix,y,i,y) =iky + y^f{x), 

d d 

^■^(a;,y;0 = + (^^ - l)/(a;), —R{x,y;^,y) =y^fix). (3.8) 

y —21 y—X 

Equations (3.7) can be simply expressed by introducing the differential operators 

^ 1 ^ ^ */i ^ d ^ Id 

These operators satisfy the same commutation relations as (2.3), i.e., [Ji, J 2 ] = [>72)^3] = *7i, 

[J3, 7i] = iJ 2 - Using (3.8) and (3.9), we can write (3.7) as 


^7(a;, y; y) = [2ikJ3 - 2f{x)Ji]T{x, y; y), 
^U(a;,y;C) = [2ikJ3 - 2f(x)Ji]L{x,y;^), 

^ A(x, y; y) = [2ikJ3 - 2f(x)Ji]R{x, y; y) + y‘^f{x). 


(3.10a) 

(3.10b) 

(3.10c) 


(We have replaced 2 ; by y.) Note that (3.10a) and (3.10b) have the same form as (2.1). 

In the same way as ( 2 . 1 ), we define the evolution operator U{x,y), which is an operator acting 
on functions of ^ and y, as the solution of the differential equation 


■^U{x,y) = [2ikJ3 - 2f{x)Ji]U{x,y) 


(3.11) 
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FIG. 3. Construction of f{x, z), L{x, z), R{x, z) from f{x, y), L{x, y), R{x, y) and r(y, z), Ri{y, z), Rr{y, z) 
[Eqs. (3.5) and (3.6)]. 



FIG. 4. Explanation of Eqs. (3.8) using diagrams. 


with the boundary condition U{y,y) = 1 (identity operator). According to general formula (2.4) 
(which is representation independent), this U can be expressed as 

U{x,y) = exp[-Rr{x,y)J+][T{x,y)]'^'^^ exp[Ri{x,y)J-], (3.12) 


where 


J+-J1+IJ2 


J- = J\— iJ2 = 




When acting on ^-independent functions of the operators J± and J 3 reduce to 


yW - _y2 A 




_a 


7(0) _ 


(3.13) 


(3.14) 


It is well known that , and 

and special conformal transformation,^*^ 


are, respectively, the generators of dilatation, translation, 
which means that 


exp(cJ 3 ^°^) 5 (C) = exp(cjl°^)g (0 = g{i + c), 










































































(3.15) 


exp(cj:[°^)5(C) = 

for an arbitrary constant c. These relations can be generalized to the /r-dependent case as 
exp(cJ 3 )g( 5 ,/x) = exp{cj-)g{^, ^i) = g{^ + c,g), 

exp{cj+)g{^,fi) = (3.16) 

The first equation of (3.16) easily follows from exp(cj73) = exp(|^^) exp(c^^). The second 
equation of (3.16) is a trivial extension of (3.15). Defining C = C a = g/^, we have 

exp{cj+)g{^,g) = exp 5(C,aC)- (3-17) 

The third equation of (3.16) follows from (3.17) and the third equation of (3.15). 

Using (3.12) and (3.16), we can see that U acts on a function g{^,fJ,) as 

l^{x,y)g{i,g) = exp(-i?rJ+)T^'^'' exp{RiJ_)g{^,g) = exp{-RrJ+)T'^'^^gi^ + Ri,fJ,) 

= exp{-RrJ+)g{T‘^^ + Ri,Tg) = g( ^ ^ | ^ ^ ^ )- (3-18) 

In view of definition (3.1), this expression reads 

l^{x,y)g{.i,g) = giL{x,y;^),f{x,y]^,g)). (3.19) 

Thus, applying the operator U{x, y) to a function g{^, g) amounts to replacing ^ and g by L{x, y; 
and T{x,y; g), respectively. In particular, 

^ix,y)g = f{x,y;^,g), U{x,y)^ = L{x,y;^). (3.20) 

When the function g in (3.19) has one or more spatial variables in addition to ^ and g, these 
variables remain unchanged on the right-hand side of (3.19). For example, 

U{x,y)g{y,^,g) = g{y,L{x,y;^),f{x,y;^,g)). (3.21) 

As can be seen from the definition, the evolution operator has the property 

U{x,y)U{y,z)=Uix,z). (3.22) 

From (3.22), it follows that U{y,x) = [U{x,y)]~^. From (3.22) and (3.20), we have 

^{x,y)f{y,z;^,g) = f{x,z;^,g), U{x,y)L{y, z; ^) = L{x, z;^). (3.23) 

This is none other than Eqs. (3.5) rewritten using U. Similarly, (3.6) can be rewritten as 

l^{x,y)R{y,Z]^,g) = Rix,z;^,g) - R{x,y]^,g). (3.24) 


IV. SOLUTION OF THE INHOMOGENEOUS EQUATION 

Unlike (3.10a) and (3.10b), Eq. (3.10c) is inhomogeneous. We define 



A=^ + 2f{x)Ju B = Jz. 

(4.1) 

and write (3.10c) as 

{A - 2ikB)R{x, y; g) = g'^fix). 

(4.2) 

More generally than (4.2), let 

us consider the inhomogeneous differential equation 



II 

1 

(4.3) 
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where h is an unknown function and g is a given function of x, and fx. 

We can rewrite (3.22) as U{z,x)U{x,y) = U{z,y). Differentiating both sides of this equation 
with respect to x, and then setting y = x and using U{x, x) = 1, we have 

x) = U{z, x)[-2ikj^ + 2/(x) J'l], (4.4) 

Therefore, for arbitrary z, 


dx 


[W(z, x)h{x, /i)] = U (z, x) 


dx 


- 2ikJ3 + 2f{x)Ji 


h{x,^,y). 


(4.5) 


Substituting (4.3), and interchanging x and z, we rewrite (4.5) as 

■^[U{x,z)h(z,^,fj,)] =U{x,z)g(z,^,fj,). (4.6) 

To solve (4.2) or (4.3), we need a boundary condition. When y is finite, the boundary condition 
for R{x,y,^, y) is R{y,y,^, y) = 0 for any ^ and y [see Eq. (3.4)]. Suppose that the solution h 
of (4.3) satisHes the boundary condition h{y,^,y) = 0, with a fixed number y, for any ^ and y. 
Integrating (4.6) from z = y to z = x gives 

Hx,^,y)=f dzU{x,z)g{z,^,y), (4.7) 

Jy 

where we have used lA{x^ y)h{y^ y) = h{y^ L(x, y\ ^), T(x, y\ y)) = 0. Hence, 

Rix,y;^,y) = f dzU{x, z)y‘^ f{z) = f dz [f{x, z; y)]^ f(z). (4.8) 

dy dy 

In this paper, we are interested in reflection coefficients for semi-infinite intervals. Corresponding 
to (1.8), the definition of R for a semi-infinite interval is 

R{x,-oo]^,y) = lim R{x,y;^,y). (4.9) 

y-^-oo 


Provided that the integral in (4.8) is convergent in the limit y —>■ —oo, we have 

R{x,-co;^,y) = f dz [f{x, z; y)]'^ f{z). (4.10) 

y —oo 

It is also possible to derive (4.10) directly by considering Eq. (4.2) with y = —oo. When dealing 
with R(x, —oo), we need to be careful about the boundary condition. The boundary condition for 
R(x, —oo) is not lima;_,._oo R(x, —oo) = 0. (This does not necessarily hold.) Setting z = —oo in 
(3.24) and using (4.9), we find that 

lim ll{x,y)R{y,-oo;^,y) =0, (4.11) 

y-s.-oo 

where x is arbitrary. This is the correct boundary condition for R{x, —oo). Let us assume that 
the solution h of equation (4.3) satisfies the same boundary condition as (4.11), 

lim U{x,z)h{z,^,y) = Q, (4-12) 

2—>• —OO 

for any ^ and y. Integrating (4.6) from z = —oo to z = x, and using (4.12), we obtain 

/ X px 

dzU{x,z)g{z,^,y)= dz g{z, L{x, z;^),f{x, z; y)). (4.13) 

-OO y —oo 


In the following analysis of R{x, —oo), we will mainly deal with functions that satisfy boundary 
condition (4.12). So, let us restrict the domain of the operator {A — 2ikB) to functions satisfying 
(4.12). Then, as shown above, the inverse of (.4 — 2ikB) is given by 

{A-2ikB)~^g{x,^,y)= f dz g{z, L{x, z;^),f{x, z; y)). (4.14) 

J —OO 

Using this inverse, the solution of (4.3) is obtained a,s h = {A — 2ikB)~^g. In other words, the 
solution of (4.3) is given by (4.14) if we assume that h satisfies (4.12). Since R{x,—oo) satisfies 
equation (4.2) and boundary condition (4.11), we can express it as 

R{x, —oo; y) = {A — 2ikB)~^y^ f{x). 


(4.15) 
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V. TRANSFORMATION OF THE VARIABLES 

The operator is diagonalized in the basis of its eigenfunctions . When dealing with the 
low-energy expansion, it is convenient (as we shall see in Sec. VII) to use a basis in which J7i, 
rather than jTa, is diagonalized. To find such a basis, we note that and are related by a 
rotation of angle 7r/2 around the 2-axis, which is represented by the operator exp(—ii7rj72). It 
follows from commutation relations (2.3) that 

exp(^yj2j Ji exp(^— —j2j=Jz, exp(^-^j2j J3exp(^-yj2j =-Ji. (5.1) 

Applying this rotation to ^ and /i gives (see Appendix B) 

5 = exp(^-y:72jC = ^ = exp(^--^:72jM = (5.2) 

Let Ji, J 2 -, Jz be defined by (3.9) with ^ and /r replaced by ^ and jx, 

~ 1 ,. X2\ ^ 1 r_ 9 ~ i ros 9 i d ~ ~ d 1,9 

If we define the variables w and a by ^ = e™ and fl = , expressions (5.3) become 

~ d \ d d % 0 

Ji = -(sinhw)---(coshw)a—, J 2 = z(coshw)--h -(sinht(;)a—, 

aw 2 aa aw 2 oa 


Thus, Jz takes a simple form. The original variables ^ and /x are related to w and a as 


^ = tanh■ 


or w = log ■ 


The first equation of (5.5) has the same form as (2.II). The variable w has the same meaning as 
in Sec. II when restricted to real numbers. (Now ^ and w are complex numbers.) 

From (5.1), it follows that Ji = Jz and Jz = —J\- Using (5.4) in (4.1) gives 




(9 1 d 

B = (sinhw)—-1—(cosh w) a—. 

aw 2 oa 


To simplify the expression for A, let us define 

W = w + V{x). (5.7) 

For any function g(x, w) = g{x, W — V(x)), we have the relations 

where {dg/dx)w and {dg/dx)w denote the partial derivatives with fixed W and fixed w, respec¬ 
tively. Let us now change the set of independent variables from {a;, w, a} to {x, IF, a}. We assume 
that the operators A and B act on functions of x, IF, and a. From (5.8), we see that the expressions 
for A and B become 

A=^, ^ = sinh[IF-F(x)]^ + icosh[IF-F(x)]a^. (5.9) 

Equations (3.21) and (4.14) can be written in the variables {x, IF, a} as 

U{x,z)g{z,W,a) = g{z,w{x, z;W), A(x, z;W,a)), (5.10) 


(A — 2ikB) ^g(x,W,a)= / (i 2 :g(z, a;(x, z; IF), A(x, 2 ;; IF, a)), 


— 00 


(5.11) 
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where 


uj{x, z] W) = log 


l + L(x,z-,W) 
l-L{x,z]W) 


V{z), 


A{x, z; W, a) 


T{x, z; W, a) 
^l-L2(x,z;W) 


(5.12) 


Expressions (5.12) come from the last two equations of (5.5) and W = w -\- V{z). [Here, W is not 
w + V{x) but w + V{z), since this W appears in g{z, W,a).] 

In (5.12), the functions T and L are written with the variables W and a. For functions which 
have two spatial variables such as T or L, the variable W is associated with the first spatial 
variable, which corresponds to the right endpoint of the interval. For example, in T(xi,X 2 \W,a), 
the definition of IF is IF = ic + V{xi). From (3.1), (5.5), and (5.7), we have 


f{x,y]W,a) 


a^{x,W)T{x,y) 
l-^{x,W) Rr{x,y)' 


R{x, y; IF, a) 


L{x,y;W) = Ri{x,y) + 


^{x,W) T'^{x,y) 

I - ^{x,W) Rrix,y)’ 


j‘^{x,W) Rrix,y) 
l-^{x,W) Rr{x,y)' 


where 

IF - V(x) , - IF - V(x) 

^(a;, IF) = tanh- - -, j{x, IF) = \/l — = sech- - -. (5-14) 

Similarly, we can rewrite Eqs. (2.12) using the variables {x,y,W} as 


T{x,y,W) 


'y{x,W)T{x,y) 

1 - ^{x,W)Rr{x,y) 

Rr{x,y;W) = 


Ri{x,y; IF) = Ri{x,y) 


^{x,W)T'^{x,y) 

1 - ^{x,W)Rr{x,y)' 


'y'^jx,W)Rr{x,y) 

1 - ^{x,W)Rr{x,y) 


ax,w), 


(5.15) 


Comparing (5.13) with (5.15), we find 


f{x,y;W,a) = aT{x,y;W), L{x,y;W) = Ri{x,y;W), 

Rix, y- IF, a) = [i?,(x, y; IF) + C(x, IF)] . (5.16) 


The variable IF is a complex number. As noted before, T, L, and R are analytic functions of ^ in 
the unit circle j^j < 1. Since F(a;) is a real number, |^(a;, IF)| < 1 corresponds to jlmlFj < 7r/2 
irrespective of x. Therefore, T, L, and R are analytic functions of IF in the strip |ImIF| < 7r/2. 


VI. EXPRESSIONS IN EIGENSPACES OF iid/dfi 


The operators 77i, 772, TTs, as well as TTl, 772, TTs, satisfy the same commutation relations as 
(2.3). So, they constitute a representation of the Lie algebra s/(2,C), where the representation 
space consists of analytic functions of w and a (or ^ and /x). This representation is reducible. All 
these operators commute with the operator defined by 


AC ^ ^ 

on aa 


( 6 . 1 ) 


When w and a are taken as independent variables, the eigenspace of M with eigenvalue u consists 
of functions of the form a''g{w). In this eigenspace, 77i and J 2 [Eqs. (5.4)] reduce to 


77 /= — sinh -coshru, 

aw 2 


= icoshw-^ 
aw 


i— sinhic, 
2 


( 6 . 2 ) 


which act on functions of w alone. That is to say, 


Ji a''g{w) = a''j}‘'^g{w) {i = I, 2, 3), 


(6.3) 
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where we have also defined = d/dw. The operators J 2 ''\ and satisfy the same 

commutation relations as (2.3). Thus, we have a different representation of sl{2, C) for each v. [In 
this representation, the Casimir operator is —-1(1 — f) times the unit operator.] 

When the spatial variable x is included and the independent variables are changed from {x, w, a} 
to {x, W, a}, the eigenspace of Af with eigenvalue i' consists of functions of the form a''g{x, W). In 
each eigenspace, the operators A and B [Eqs. (5.9)] act as 

Aa''gix,W) = a''A^''^gix,W), Ba''g{x,W) = a''B^‘'^g{x,W), (6.4) 

A^'''> = —, B^'''> = sinh[IT - V(x)]^ + - cosh[IT - V(x)]. (6.5) 

OX aW 2 

The inverses of {A — 2ikB) and — 2ikB^^^) satisfy the corresponding relation 

{A - 2ikB)-^a''g{x, W) = a''{A^''^ - 2ikB^''^)-^g{x, W). (6.6) 

Using (5.II) for the left-hand side of ( 6 . 6 ), and inserting (5.16), we find 

{A^'''^-2ikB^''^)-'^g{x,W)= dz -—— g{z,ui{x, z;W)). (6.7) 

J_oo [I- Rt{x,z;W)\ 

If {xA,b} are taken as the independent variables instead of {x, IT, a}, the eigenspace of Af 
consists of functions of the form ijyg{x,^). From (4.1) and (3.9), we have 

■Ag!"g{xA) = ^ a'" = a’"^(v)9{xA), ( 6 - 8 ) 

As shown in (5.5), g = 'ya with 7 = \/l—Comparing (6.4) with ( 6 . 8 ), we can see that 

7 ", %) = 7“" (6.10) 

The number iz corresponds to the number of white circles in each diagram in Fig. 2. The 
functions T, L, and R belong to the eigenspaces with = 1, 0, and 2, respectively. 


VII. LOW-ENERGY EXPANSION 


Now let us study the expansion of (4.15) in powers of k. Here, we assume that the following 
limits exist: U(— 00 ) = lim 2 ;_>_oo U(x), /(— 00 ) = lim 2 ,_,._oo/(a:), and /'(— 00 ) = lima;_ 5 ._oo/^(a;). 
(In the last expression, /' means the derivative of /.) These limits may be infinite. 

For A: = 0, we can easily solve Fq. (1.5) and obtain 


x{x, y; k = 0) = cosh 


V{y) -V{x) 


l3{x, y; k = 0) = sinh 


V{y) - V{x) 


Substituting (7.1) into (1.7), (5.13), and (5.12) successively, we find 

U!(x, y; IF, a; fc = 0) = IF, A{x, y; W, a; k = 0) = a. 
Setting fc = 0 in (5.11), and inserting (7.2), we obtain 


A ^g{x,W,a)= ( g{z,W,a)dz. 

J —00 


Similarly, on setting fc = 0, condition (4.12) becomes lima;^_oo h{x, IF, a) = 0 [see also Fq. (5.10)]. 
Indeed, the inverse of 7l = d/dx is given by (7.3) if we assume that the domain of A consists of 
functions that vanish as a; —>■ — 00 . 
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Let g be a function for which A~^g makes sense. Since AA~^g = g, we have the identity 

g = {A — 2ikB)A~^g + 2ikBA~^g. (7-4) 

The inverse of {A — 2ikB) is given by (5.11) if the domain of {A — 2ikB) is restricted to functions 
satisfying (4.12). This means that {A — 2ikB)~^{A — 2ikB)A~^g = A~^g if 

lim I4{x, z){A~^g){z) = 0. (7.5) 

z—^ — oo 

Let us apply {A — 2ikB)~^ to both sides of (7.4). If (7.5) holds, we have 

(A - 2ikB)-^g = A-'^g + 2ik{A - 2ikB)-^BA~'^g. (7.6) 

The expressions {A — 2ikB)~^g and A~"^g make sense if and only if the integrals in (5.11) and (7.3) 
are convergent. The second term on the right-hand side of (7.6) automatically makes sense if both 
{A — 2ikB)~^g and A~^g make sense. In summary, (7.6) holds if 

(i) the integrals in (5.11) and (7.3) are both convergent and 

(ii) condition (7.5) holds. (7.7) 

The -I- 1 times iteration of (7.6) yields 

N 

(A - 2ikB)-^g = '^{2ikT{A~^BTA-^g + {2ik)^+^{A - 2ikBy^B{A-^B)^A~^g. (7.8) 

n—0 

We assume that {A — 2ikB)~^g makes sense. Then (7.8) holds if {A~^B)'^A~^g makes sense for 
all n < IV and if 

lim U{x,z)[{A~^B)'^A~^g]{z) = Q for all n < A^. (7.9) 

2—>• —OO 

From (7.8) and (4.15), we obtain 

R{x, —oo) = ro -I- ikri + (ik)'^r 2 + {ik^r^ H--I- {ik)^fN + Pn (7-10) 

with 

rn{x) = {2A~^B)'^A~^pL'^f{x), (7.11) 


Pn{x, k) = 2{ik)^^^{A — 2ikB)~^BfN{x). 

(We omit to write the dependence on W and a.) From (7.9), we have the condition 

lim lA{x,z)fn{z) = 0. 

2—>• —OO 

Equation (7.10) is valid if 

(i) the right-hand side of (7.11) makes sense for n = 0,1,..., A^ and 

(ii) condition (7.13) holds for n = 0,1,..., A^. 

Substituting p = aj cosh using (6.4), we can write (7.11) as 

f„(i) [2(.4W)--BW]"(.4g>)-' 

where is given by (6.5) with ly = 2, and 

iA^^Y^g{x,W)= r g{z,W)dz. 

J —OO 

Since f{x)/ cosh^ w-V(x) _ tanh obtain from (7.15) as 


fo(x) = 


W-V(x) 


— tanh 


W-V{-oo) 


(7.12) 

(7.13) 

(7.14) 

(7.15) 

(7.16) 

(7.17) 


We need to consider three cases, E(—oo) ±oo, E(—oo) = -boo, and 1L(—oo) = —oo. 
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A. The case V(—oo) = Vi ^ ±oo 

First, we consider the case where 14(— 00 ) has a finite value li. We can rewrite (7.17) in the 
form 


ro(x) = 




cosh cosh 


(7.18a) 


Using r„ = [ 2 (yl( 2 ))-i^( 2 )]" ‘ fo, the coefficients for n > 1 can be calculated as 


31 


fi(a;) = 


r2(x) = 


a 


cosh 


2 W-Vi 


dz sinh[li — V(z)], 


a 


cosh 


3 W-Vi 


L{w+v,)/ 2 f j rf2'e-^(^)sinh[Ui-U(z')] 

I J —QO J —00 

+ [ dz f dz' sinh[Ui - V(z')] \, 

J —00 J —00 I 


and so on. It can be easily shown by induction that r„ (n > 1) has the form 


n—1 


fn{x) = H exp[(j - IF] Cj(x), 


2 j=0 


(7.18b) 


(7.18c) 


(7.19) 


where Cj{x) is a function that vanishes as x —>■ — 00 , provided that r„ makes sense. Let F, 
denote the set of real-valued functions g{x) such that 


(-) 


/ XQ 

(1 -b |x|")|( 7 (x)| dx < 00 for any finite xq. 

-00 


(7.20) 


It is not difficult to show that r„ makes sense if U — Ui G F^_\ (see Appendix C of Ref. 26). Since 
C F^Z\ for n < N, condition (i) of (7.14) is satisfied if U — li G 
We also need to check condition (ii) of (7.14). From (5.10), (5.12), and (7.18a) we obtain 


U{x,z) fo{z) 


—r^(x, z) tanh 
1 + L{x, z) tanh 


(7.21) 


Similarly, from (5.10), (5.12), and (7.19) we have, for n > 1, 


U{x, z) fn{z) 


f\x,z) 


cosh 


n+1 V{z 


— ^exp[(j 


1=0 


^)V{z)\B,{x,z)C,{z), 


(7.22) 


where 


Bj(x,z)= 1 + L{x,z)Y \l — L{x, z) ” ^ 


1 -b L{x, z) tanh 


-(n+l) 


(7.23) 


The asymptotic behavior of T{x,z) and L{x,z) as z —1 —00 is shown in Appendix C. As we are 
assuming that /(— 00 ) exists, U(— 00 ) ^ ±00 implies /(— 00 ) = 0. From (C6) of Appendix C, 
we can see that limz_>_oo \Bj{x,z)\ < 00 . Since tanh q Cj{z) —>■ 0 as z —>■ — 00 , 

the right-hand sides of (7.21) and (7.22) vanish in this limit. This means that (7.13) holds as 
long as fn makes sense. Both conditions (i) and (ii) of (7.14) hold, and hence, (7.10) is valid, if 
u - Ui G 
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B. The case V(—oo) = +oo 

Next, we study the case ^(— 00 ) = + 00 . Now (7.17) becomes 

ro(x) = 2 a 2 + 

Using r„ = [ 2 (^( 2 ))-i^( 2 )j" ro, we can calculate 

fi(x) = 20^6^ r dze-^^^\ hix) = T dz T dz' 

J —00 J —00 J —00 

and so on. For general n > 1, the expression for r„ can be explicitly obtained as^® 


i(x) = 


E E ■ 

(72—^1 (Ts — H 


■ E 

<7 n — i 1 


exp 


(i + E'^j) 

i=2 


W 


f+l,tT2,...,cr^ 


.(a;), 


where 



(7.24a) 


(7.24b) 


(7.25) 


(7.26) 


n J 

C(Tl = 2 1 ^^^^ (Ti^ CTj . 

J=1 i=l 

It can be seen from (7.27) that C+i^cy 2 ,...,<T„ = 0 unless 

j 

1 + ^ CTi > 0 for j = 2 ,3,..., n. 

i=2 


(7.27) 


(7.28) 


So, only the terms with a 2 ,---,cr„ satisfying (7.28) are present in (7.25). [In particular, only 
(72 = +I survives in (7.25).] It can be shown^® that |/+i,cr 2 ,...,<T„ | < 00 for all CT 2 , ..., Cn satisfying 
(7.28) if e~^ G Pi-\- Therefore, rn makes sense if e“^ G and condition (i) of (7.14) is 

satisfied if e~^ G 

We proceed to check condition (ii) of (7.14). Using (5.10) with (7.24a) and (7.25), we obtain 

T‘^ (X z') 

U{x,z)r,{z)= r (7.29) 

1 J-J y ^ / 


U{x,z)fn{z) =f‘^{x,z) ^ ^ 

(T2=±1 (T3 = ±1 


E ^+1 ,(T2,...,<7n F(J2 ,...,(Tn {X^ Z') 77^2 . .,(T„ ('^) J 

cr„=±l 


(7.30) 


where 


B(j 2 ,...,<Tn {x, z) — 1 + L(x, z) 


1:7=2 r 


Br 


n 

= exp (l + ^crj)u(z) 


1 — L(x, z) 

7+1,(T2,...,o-„ (z)- 


2-E7=2 


J =2 


(7.31) 

(7.32) 


When U(— 00 ) = + 00 , there are three possible cases for /(— 00 ), namely, /(— 00 ) = + 00 , 
/(— 00 ) = c > 0, and /(— 00 ) = 0. Let us first assume that /(— 00 ) ^ 0. From (Cl) and (C3) of 
Appendix C, we know that if /(— 00 ) = +00 or /(— 00 ) = c > 0, 

T(x, z) —>■ 0, L{x, z) —>■ constant ^ 1 as z —>■ — 00 . (7.33) 

(Here, we are assuming that c > jfcj, as we are now interested in the low-energy region.) In this case, 
we can easily show that limz_>_oo |71cr2,...,<T„ {z)\ <00 for all CT 2 , ..., cr„ satisfying (7.28). On account 
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of (7.28), we have Yl^= 2^3 — 0- From this and (7.33), it follows that lim^^-oo \Bcr 2 ,...,a„{x, z)\ < 
oo. Now it is obvious from (7.33) that the right-hand sides of (7.29) and (7.30) vanish in the limit 
z —>■ —oo. Thus, (7.13) holds if f„ makes sense. 

The situation is different for the case /(—oo) = 0. When t^(—oo) = -l-oo and /(—oo) = 0, 
condition (7.13) does not hold for Imk = 0. As can be seen from (C4), in this case, T{x,z) 
oscillates as z —>■ —oo for Im/c = 0, so neither (7.29) nor (7.30) becomes zero in the limit z —>■ —oo. 
However, this does not mean that (7.10) is not correct for Imfc = 0. Corresponding to (1.9), the 
definition of R for real k is R{k) = limej,o R{k+ie). If conditions (7.14) hold for Im/c > 0, Eq. (7.10) 
is valid even for Imfc = 0 as long as remainder term (7.12) is interpreted as limg^o Pn{x, k + ie). 

Let us assume Imfc > 0 to consider the case where y(—oo) = -l-oo and /(—oo) = 0. It can be 
seen from (C4) that T{x,z) approaches zero exponentially as z —>■ —oo. Meanwhile, Da- 2 ,...,(t„{z) 
defined by (7.32) becomes infinite in this limit; it grows at most like |z|” as z —>■ —oo if makes 
sense. Since T^{x,z) falls off exponentially and \Ba 2 ,...,a„{x, z)\ remains finite, both (7.29) and 
(7.30) vanish in the limit z —>■ —oo. 

Thus, in the case y(—oo) = -l-oo, Eq. (7.13) holds (at least for Imfc > 0) if e~^ G 

irrespective of whether /(—oo) 0 or /(—oo) = 0. If e~^ G both conditions (i) and (ii) of 

(7.14) hold for Imfc > 0, and hence, (7.10) is valid for Imfc > 0. 


C. The case M(—oo) = —oo 

The case E(—oo) = —oo is similar to the case V{—oo) = -l-oo. In this case we have 


= H ^ C'-i.<^ 2 .....^n exp 

(T2 = il (T3 = ±1 


n 

(-1 + Y,<^3)w]l-l,a,,...,aAx). (7.34) 


i=2 


In the same way as in the previous case, we can show that (7.10) is valid if e'^ G 


Thus, we have found that condition (ii) of (7.14) is always satisfied as long as condition (i) holds. 
Equation (7.10) is valid if the coefficients fg, fi,..., tat make sense. 

For (7.10) to be an asymptotic expansion, the remainder term must satisfy 

lim/5Ar(a;, A:)/A:'^ = 0. (7.35) 

k^O 

Using (5.11), the expression for remainder term (7.12) can be written as 

PN = 2{2ik)^+^ f dz ^ KNiz,uj{x,z;W)), (7.36) 

j —OO 1 — (x, z] W J 

where Kn{x,W) = {1 /B f n. We can investigate the validity of (7.35) by using this integral 
form. This was essentially done in Refs. 24 and 26. It can be shown that (7.35) is satisfied if 
U — Ui G or G (See Appendix E of Ref. 26 for details.) 


VIII. ASYMPTOTICALLY PERIODIC POTENTIALS 

In Sec. VII, it was assumed that U(—oo) exists. An important class of potentials that does 
not fall into this category is asymptotically periodic potentials,i.e., potentials that tend to a 
periodic function as x —>■ —oo. [Here, the term “potential” refers to the Fokker-Planck potential 
V {x). The corresponding Schrodinger potential Vs (x) is also asymptotically periodic under the 
following assumptions.] In this section, we consider such potentials. We assume that V and / have 
the forms 


^(x) = Fp(x) -k Va(x), f(x) = fp(x) + fA(x), 


(8.1) 
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where Vp and /p are periodic functions with period L, while Va and /a vanish as a; —>■ —oo, 

Vp{x + L) = Vp{x), fp{x + L) = fp{x), lim Va(x)= 0, lira fA{x) = 0. (8.2) 

x—y — oo x—¥ — oo 


We also assume that the derivative of fAix) vanishes as x —>■ —oo. 

For such V{x), expression (7.17) does not make sense. This is because (7.3) does not give the 
inverse of A when the domain of A consists of functions that do not vanish as x —>■ —oo. We need 
to find the appropriate inverse of A for asymptotically periodic functions. 

Let us consider Eq. (4.3) with k = 0, 


Ah{x, W, a) 


dx 


h{x, W, a) = g{x, W, a), 


assuming that both h and g are asymptotically periodic functions of x, 


h — hp + /^Aj 


5 — 5p + 5 Aj 


(8.3) 


(8.4a) 


hp{x + L) = hp{x), gp{x + L) = gp{x), lim hA{x)=0, lim gA(a;) = 0. (8.4b) 

x—¥ — oo X—y — OO 

We also assume that h and g are analytic functions of W in the region |ImlF| < 7r/2. (See the 
comment at the end of Sec. V.) Equations (8.3) and (8.4) imply gA = dhAjdx and gp = dhp/dx. 
Hence, it follows that 


c—L 


gp{z,W,a) dz = 0 


for any x. The relation gA = dhA/dx can be inverted as 

r 

hA{x,W,a)= / gA{z,W,a)dz. 

J —OO 


(8.5) 


( 8 . 6 ) 


On the other hand, gp = dhp/dx cannot be uniquely inverted. All we have is 

hp{x,W,a)= f gp{z,W,a)dz + Cxo(W,a), (8.7) 

Jxo 

where xq is an arbitrarily fixed number, and CxoiW,a) is an x-independent function which is 
undetermined at this moment. 

When g is given, the function h satisfying (8.3) is not uniquely determined by conditions (8.4), 
as Cxp{W,a) is undetermined. However, assuming that we have a further condition to determine 
Cxo{W,a) (such a condition will be shortly shown), the inverse of A can be written as A~^g = 
hp + Ha- We denote the right-hand side of (8.7) by Ap^gp{x,W,a). Then, 


A ^g{x,W,a) 


9 A{z,W,a) dz- 


Ap^gp{x,W,a) 


( 8 . 8 ) 


It should be noted that A~^g does not make sense unless gp satisfies (8.5). 

Using this A~^, we shall derive the expansion of {A—2ikB)~^ as we did in Sec. VH. The starting 
point is Eq. (7.6). Although Cx^iW, a) in (8.7) has not yet been determined, Eq. (7.6) holds for any 
Cxo as long as condition (7.5) is satisfied. However, for (7.6) to be useful for low energy expansion, 
it is necessary that 


\\-m 2ik{A — 2ikB) ^BA ^<7 = 0 (8.9) 

fc—s-o 


so that (7.6) may yield limfc_j.o(A — 2ikB) ^g = A ^g. We also require that 


YiuY{A — 2ikB) ^BA ^g = A ^BA ^g. 


( 8 . 10 ) 
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Then (8.9) holds if A~^BA~^g makes sense. In order that A~^BA~^g may make sense, the periodic 
part of BA~^g must satisfy the condition corresponding to (8.5). The periodic part of BA^'^g is 
B-pA~^g-p with B-p defined by 

Bp = sinh[lT - + 2 cosh[lT - Vp(x)] a—. (8.11) 

Putting BpAp^gp in place of gp in (8.5), we have 


/ x — L 


{BpA '^gp){z, W, a) dz = 0. 


( 8 . 12 ) 


The function Cxo(W,a) is determined by condition (8.12). 

It is easier to work in the eigenspace of A/” and deal with A^‘''^ and defined by (6.5). We 
assume that g{x, W) = gp{x, W) + g\{x, W) as in (8.4). Corresponding to (8.8), we have 

{A^''^)-^gix, W)= ( gA(z, W) dz + (A^i^'^y^gpix, W), (8.13) 

J —00 

{A^;y)-^gp{x,W)= r gp{z,W)dz + C^:J{W). (8.14) 

J Xq 

The function C'!ig iW) in (8.14) is determined from the condition corresponding to (8.12), 

r [B^^\AP)-^gp\{z,W)dz = 0, (8.15) 

J X — L 

where 

d ly 

B^p'^ = sinh[lT - ^ cosh[lT - Tp(a;)]. (8.16) 

Note that {A^''^)~^ explicitly depends on although the expression of A^''^ in (6.5) appears to be 

independent. This is because the appropriate domain of A^'''^ depends on v. 

For the analysis of R, we need only to deal with v = 2. Let us determine {W) from condition 
(8.15). For V = 2, substituting (8.14) and (8.16) into (8.15) gives 


We define the constants Lq and Vb by 
Lq = 


dz'gp{z',W) + C^^Kw) 


= 0 . 


Co ^ t log 


1 , 


j:_,e-yA^)dz^ 


so that e^yp^^'>dz = Lq 6^^°. Integrating (8.17) with respect to W, we obtain 




-1 


Lq sinh(lF — Vg) 


(8.17) 


(8.18) 


f dz f dz'sinh[LF — Cp(z)](/p(z', LF) + cl , (8.19) 

J X—L 'J xq ) 


where C is a constant of integration. By our assumption, Cxq yW) should be analytic in the strip 
|ImlF| < 7r/2, so the right-hand side of (8.19) should not be singular at IF = Vq. From this 
condition, the constant C in (8.19) is determined as 


C = - 




Fp(z)]ffp(z',Fo). 


( 8 . 20 ) 


Since xg is arbitrary, we may let xg = x — L. Substituting (8.19) into (8.14), letting xg = x — L, 
and using gp{z) dz = 0, we obtain 




f dz f dz'{sinh[Fp(z) - lF]gp(z', IF) - sinh[Fp(z) - Fo]gp(z',Fo)}. (8.21) 

JX—L Jx—L 
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Thus, with (-4p^^) ^ given by (8.21), the inverse of cam be written as 

(A^^Y"g{x,W) = {A^^Y"gp{x,W)+ [ g^{z,W)dz. (8.22) 

J —OO 


The same expression as (8.21) was derived in Ref. 2 by a different method. Although the method 
of Ref. 2 is meaningful for its own sake and for its applications, expression (8.21) itself can be 
derived more simply from condition (8.15) as we have seen above. 

We can calculate the coefficients f„ of (7.10) by using (8.22), instead of (7.16), in (7.15). The 
calculation is the same as that for the expansion of Rr studied in Ref. 27. We obtain 


ro(x) = tanh 


fi(x) = 


W-V(x) 


— tanh 


W-Vo 


■ sinh 


Vo-Vix) 


coshl^ coshi^^^M’ 


cosh^ 


/ dz'smh[Vp{z')-Vp{z)] 

y-^-^0 Jx—L Jx — L 

+ J dz {sinh[Vb - V(z)] - sinh[Vb - Vp( 2 ;)j}'j , 


(8.23a) 


(8.23b) 


and so on. It can be shown^^ that f„ makes sense if Va G Fi-\- So, condition (i) of (7.14) is 
satisfied if Va G ^atIV When f„ (n > 1) makes sense, it has the form 


r„(x) = ^ H,{W) [q{x) + Cf (x)] , (8.24) 

i=i 


where m is a finite number, each C']’(x) is a periodic function, each C^{x) is a function that 
vanishes as x —>■ —oo, and each HjiW) is an analytic function of W in the strip |ImIT| < 7r/2. 
[Note that (8.23b) has the form of (8.24) with m = 1.] 

Let us check condition (ii) of (7.14). From (8.23a), (8.24), and (5.10), we obtain 


and, for n > 1, 


U{x,z) ro{z) 


—r^(x, z) tanh 
1 + L{x, z) tanh 


(8.25) 


U(x,z) fn(z) 


f^(x,z;W,a)J2 

j=i 


Hj(uj(x,z;W)) 

l-P(x,z;W) 


[C^{z) + Cf{z)]. 


(8.26) 


Unlike (7.21) and (7.22), the right-hand sides of (8.25) and (8.26) do not tend to zero as z —>■ —oo 
if Im/c = 0. In this case, we need to assume Ivcik > 0. As mentioned in Sec. VII, this is enough 
for the validity of (7.10) in the closed half plane Imfc > 0. 

Let Up(x, y; k) be the 2x2 matrix defined as the solution of (1.5) with / replaced by /p, satisfying 
the boundary condition Up{y, y; k) = I. Let A be the eigenvalue of f7p(x, x —L; k) such that |A| > 1. 
(Such A uniquely exists for Im fc > 0, since the two eigenvalues of Up are reciprocal of each other. 
This A is independent of x. See Ref. 2 for details.) The asymptotic form of T for Im/c > 0 can be 
written in terms of this eigenvalue A, 

T{x,z]W,a) ^ aY^^C{x,z-,W) as z —>■ —oo, (8.27) 


where C(x,z]W) is a periodic function of z. This implies that T{x,z) —>■ 0 as z —>■ —oo, since 
|A| > 1. The asymptotic form of L is 


L(x, z; W) ~ Rf{oo, z) 


I3p{z + L,z-,-k) 
A~i — ap(z -I- L, z; k) 


as z —>■ —oo. 


(8.28) 


where Op and /3p are the elements of the matrix Up defined in the same way as in (1.6), and 
Rf = — j3p{—k) / ap{k) as in (1.7). [See Eq. (D.7) of Ref. 2.] This Rf{oo,z) is a periodic function 
of z, and |i?f(oo,z)| < 1 for Imfc > 0. 
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From (8.28) and (5.12), we can see that both L{x, z) and uj{x, z; W) tend to periodic functions 
as z approaches —oo. Hence, each term in the sum on the right-hand side of (8.26) tends to a 
periodic function. These terms are always finite. Since T{x,z) —>■ 0 as z —>■ —oo, the right-hand 
sides of (8.25) and (8.26) vanish in this limit. Thus, (7.13) holds if makes sense. In this case, 
too, condition (ii) of (7.14) is satisfied if condition (i) holds. 

It can also be shown^”^ that (7.35) holds if Va G Therefore, for asymptotically periodic 

potentials, (7.10) is valid as an asymptotic expansion to order if Va G 

IX. HIGH-ENERGY EXPANSION 

In this section, we study the expansion of (4.15) for large |fc|. As |fc| becomes large, it is expected 
that {A — 2ikB)~^ will approach —{2ik)~^B~^ in some sense. However, the inverse of B does not 
exist if we let the domain of B include all analytic functions of ^ and /i. Instead, we confine 
ourselves to the eigenspace of Af with eigenvalue > 0, in which B reduces to defined by 
(6.9). We let the domain of consist of functions of ^ which are analytic in the region |^| < 1. 
[See the comment above Eq. (3.3).] Then, the inverse of B(^^) exists as shown below. Suppose that 
F(u)h = g. With (6.9), this equation reads 

+ 0 HO = = 5(C)- (9.1) 

Hence, h{0 = where is a constant. As h{0 is analytic in j^j < 1, 

this ^0 must be zero. [Otherwise, h{0 becomes singular at ^ = O.j Thus, the inverse of B(^^) is 

S(~) 5(C) = 0 C-'+^"/"^5(C) dO (9.2) 

Here, we take {x, 0 b} as the independent variables rather than {x, W, a}. Eigenfunctions of Af 
have the form g''g{x, 0- Since B 5 = g = (Ag, we can write 

H 9{x,0 =-^[{■^-‘^ikB)-A]0'B'^J^g{x,0- (9.3) 

We have {A — 2ikB)~^{A — 2ikB)gyB^^^^ g = 5 if 

lim U{x,z)iA{B7\g){z,O = 0- (9.4) 

2—>• —OO 

We assume that {A — 2ikB)~^g,'"g makes sense. Letting {A — 2ikB)~^ act on (9.3), we obtain 

(A - 2ikB)-HH = 5 + ^{A - 2ikBrH''A^^)B0] g (9.5) 

provided that (9.4) holds. [We have also used (6.8).] The N times iteration of (9.5) gives 

N-l 

(A - 2.kBr7-9 = ->‘‘'0 i2-tpTr(eM-^(.))”8M » 

+ (Sir ^ 9. (f-O 

Equation (9.6) holds if B'^J^ g makes sense for all n < TV — 1 and if 

(z,C) = 0 for all n < TV - 1. (9.7) 

Using (9.6) in (4.15) yields 

-"> = 2S'‘ + (sW® (W'’ ^ (Sir ” 



with 
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c«(a;) =-/x^(S(2)^(2))" V(a;), (9.9) 

SNix,k) = j^-^{A-2ikB)~^AcN- (9.10) 

[We used f{x) = f{x) in (9.9).] From (9.7), we have the condition 

lim ll{x,z)cn{z) = 0. (9.11) 

2^ — 00 

Expression (9.8) is valid if 

(i) the right-hand side of (9.9) makes sense for n = 1, 2,..., and 

(ii) condition (9.11) holds for n = 1, 2,..., iV. (9.12) 

From (6.9) and (9.2), we have 

= j ^9{xA)d^ + ^^[{^-^^)9{x,0-9{x,Q)\- (9.13) 

Using (9.13), we can calculate (9.9) as^^ 

Cl = C2 = (-/' + £3 = (-/" + f + 2//'C - fe)9\ 

54= [-/'" +5/V'+ (/'" +2//"-2/4)5-3/Ye"+ /"eV, etc. (9.14) 

(The primes denote derivatives with respect to x.) In general, c„ has the form 

n—1 

c„ = /i2^FY4, (9.15) 

3=0 


where each Fj is a polynomial in /, /', /",..., It is obvious that c„ makes sense if / G C”“4. 

Therefore, condition (i) of (9.12) is satisfied if / G . 

We need to examine condition (ii) of (9.12). Equations (9.15) and (3.21) give 

n—1 

U{x,z)cn{z) = [f{x,z)]‘^'^Fj{z)[L{x,z)y. (9.16) 

3=0 

By using (9.16) and (C1)-(C4) of Appendix C, we can easily check whether (9.11) holds or not. 
When /(—oo) exists and is finite, (9.11) holds for Imfc > 0 if 

lim e“®/^’”^(a;) = 0 for any a > 0 and for any 0 < m < n — 1. (9.17) 

x—¥ — oo 

(Here, denotes the mth derivative of /.) Equation (9.11) holds for Imfc > 0 if (9.17) is 
satisfied for a = 0 instead of a > 0. When /(—oo) = ±oo, (9.11) holds for Imfc > 0 if 

lim = 0 for any a > 0 and for any 0 < m < n — 1. (9.18) 

X—¥ — 00 

When the potential is asymptotically periodic as in Sec. VIII, Eq. (9.11) holds for Imfc > 0 if 
(9.17) is satisfied. This can be seen by using (9.16), (8.27), and (8.28). 

The two conditions of (9.12) are satisfied for Imfc > 0 if / G and if either (9.17) or (9.18) 

holds with n = N. Then (9.8) is valid for Imfc > 0. [For Imfc = 0, remainder term (9.10) should 
be interpreted as limej,o S{x,k + *e)-] 

For (9.8) to be an asymptotic expansion, it is also necessary that 

lim k^SNix,k) = 0. 

\k\—¥(yo 


(9.19) 
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Using (4.14), remainder term (9.10) can be expressed as 


^N{x,k) = J dz[f{x,Z]k)fJ2Gj{z)[L{x,z;k)y, (9.20) 


where each Gj is a polynomial in /, /', /"..., For fixed x and z, we have 

f(a;,z;fc) i(x, z; A:) ~ g2*fc(^-z) as |A:| ^ oo. 


(9.21) 


Using (9.21), and also using the dominated convergence theorem and the Riemann-Lebesgue 
lemma, we can show that (9.19) holds in the half plane Imfc > e with arbitrary e > 0 if^^ 


/ X 

e^^\y^\z)\dz < oo for any a > 0 and for any 0 < m < N — 1. 

-CO 


(9.22) 


In particular, (9.19) holds in the half plane Imfc > 0 including the real axis if (9.22) holds for 
a = 0 instead of a > 0. Condition (9.22) is not satisfied if f{x) diverges exponentially or faster 
as X —>■ —oo. In such cases, (9.19) still holds in the sector e < argfc < tt — e with e > 0 if the 
derivatives of f{x) are sufficiently well behaved as x —>■ — oo. 


X. CONCLUDING REMARKS 

In this paper, we studied low-energy expansion and high-energy expansion of the reflection 
coefficient Rr(x, —oo) by considering its generalized form i?(x, —oo;/r, ^). The introduction of the 
two additional variables [i and ^ gives a natural description of transmission and reflection processes 
and enables us to see clearly the meaning of the formulas used for the expansions. In particular, the 
most important formula is Eq.(5.11) [or, equivalently, (4.14)]. This basic formula can be derived 
quite naturally in our generalized framework. 

The low- and high-energy expansions of Rr{x, —oo) are obtained from (7.10) and (9.8) by setting 
/X = 1 and ^ = 0 [or a = 1 and W = V{x)]. These asymptotic expansions are justified if (i) 
each coefficient of the expansion makes sense, (ii) each coefficient of the expansion belongs to the 
appropriate domain of the operator A — 2ikB, and (iii) the remainder term is of higher order in the 
expansion. We have rigorously examined these conditions. Condition (ii), especially, is essential 
to the validity of the expansions. (Conditions (i) and (iii) have been discussed in the previous 
works.^^“^^) For the low-energy expansion, we have found that in all cases under consideration, 
condition (ii) holds if condition (i) is satisfied. For the high-energy expansion, condition (ii) is not 
guaranteed by condition (i). 

The formalism introduced in this paper can be extended to a more general description of one¬ 
dimensional problems. The operator Li defined by (3.11) belongs to an infinite-dimensional repre¬ 
sentation of the Lie group SL{2, C), with a representation space consisting of functions of /i and 
We can extend this representation to SL{3, C). By utilizing this SL{3, C) structure, it is possible 
to derive the expansions of the Green function directly, without using the reflection coefficients. 
This method will be exploited elsewhere. 


Appendix A: DERIVATION OF (3.8) 

Let e = x — yhe small. To first order in e, the diagrams contributing to T{x, y; y) are the ones 
shown in Fig. 4(a). In the second diagram, the integration over the position of / is implied. The 
contribution from these two diagrams is 

r 

y explike) + / f{z)exp[ike + 2ik{x — z)]dz = y + ikye + y^f{x)e + o{e). (Al) 

dy 

Thus, r(x, y;^, y) = y + [iky -I- /i^/(x)]e -f o(e). Since T(x, x; y) = y, we obtain 


■^T{x,yA,y) 


y^x 


= lim 

£->o e 


(A2) 
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This is the first equation of ( 3 . 8 ). Similarly, the remaining two equations are obtained from 

Lix, y;^) =^ + [2ik^ + - l)/(a:)]e + o(e), R{x, y; /x) = y‘^f{x)e + o(e), (A 3 ) 

which can be seen from Fig. 4 (b) and Fig. 4 (c). 


Appendix B: DERIVATION OF (5.2) 


Let us define ^0 = (^cos| — sin|)/(cos| +^sin|) and ycg = /i/(cos| +^sin|). It is easy to 
show that 



2 




d 1 + a Cm ^ rr 

=- 


(Bl) 


Since Ce=o = C M 0 =o = M: from (Bl) we have exp(f 6 *i 72 )C = Ce exp(z 0 j 72 )/i = ye- Setting 
6 = - 7 r /2 yields exp(^ J2)C = (1 + 0/(1 “ 0 and exp(^ J2 )m = V^m /(1 “ 0 - 


Appendix C; ASYMPTOTIC BEHAVIOR OF f{x,z) AND L{x,z) AS z ^ -00 


The behavior of T{x,z) and L{x,z) as z —>■ —00 differs according to the behavior of f{x) as 
X —>■ —00. As in Sec. VII, we assume that /(—00) and /'(—00) exist. We also assume that 
\imx^-oo f (x)/ f{x) exists. (This limit may be infinite.) We have the following asymptotic ex¬ 
pressions of T and L as z —>■ —00. (For details, see Appendix F of Ref. 26 .) 

If /(—00) = 4-00, 


T{x, z) = a C\ exp 
If /(-oo) = -00, 

T(a;, z) = aC2 exp 
If /(-oo) = c ^ 0 , 

T{x, z) = a C3 exp 


“2^(0 +’ll (^) 


^V{z) + y2{z) 


[1 + 0(1)], L{x,z) = -l + o{l). (Cl) 


[l + o(l)], T(x,z) = 1 + 0(1). (C 2 ) 


a/c^ — z + ir]3{z)\ [1 + 0(1)], L{x, z) = ^ +o(l). (C 3 ) 


If /(-oo) = 0 , 


T{x, z) = aCi exp[—ffcz + f?74(z)] [1 + o(l)], L{x, z) = C5 exp[— 2 zA:z + 2z?74(z)] + o(l). (C 4 ) 


In the above expressions, each Q (z = 1 , 2 , 3 , 4 , 5 ) is a function of x, IF, and k which does not 
depend on z, and each rji (z = 1, 2, 3 , 4 ) is a function of z and k such that 


rii{z) = o(|z|) as z —>■ —00. 


(C 5 ) 


When Im/c = 0 , the function r]4{z) in (C 4 ) is real valued, and jCsj < 1 . So it follows from (C 4 ) 
that if /(—00) = 0, 


lim |T'(x,z)| = 0 forlmfc> 0 , 

—OO 

lim |L(x,z)| = 0 forlmfc>0. 


z—^ — oo 


lim [/"(x, z)| = [0(741 < 00 forlmfc = 0, 

■ —^ — 00 

lim |L(x, z)| = |(75| < 1 forlmfc = 0 . (C6) 


z—> — oc 
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